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Abstract 

Neutrino masses and mixings are analyzed in terms of left-handed 
fields and a 6 x 6 complex symmetric mass matrix M. whose 
singular values are the neutrino masses. An angle 9 U characterizes 
the kind of the neutrinos, with 9 V = for Dirac neutrinos and 
6 U = 7r/2 for Majorana neutrinos. At 9 V = baryon-minus-lepton 
number is conserved. If 8 U w 0, the six neutrino masses coalesce 
into three nearly degenerate pairs. Thus the tiny mass differences 
exhibited in the solar and atmospheric neutrino experiments are 
naturally explained by the approximate conservation of B — L. 
Neutrinos are nearly Dirac fermions. 

This B — L model leads to these predictions: neutrinos oscil- 
late mainly between flavor eigenfields and sterile eigenfields, and 
so the appearance of neutrinos and antineutrinos is suppressed; 
neutrinos may well be of cosmological importance; in principle 
the disappearance of v T should be observable; and Qu(3(3 decay is 
suppressed by an extra factor of 10 -5 and so will not be seen in the 
Heidelberg/Moscow, IGEX, GENIUS, or CUORE experiments. 
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The current wisdom on neutrinos is that the seesaw mechanism [1] forces their 
masses to be very small. This paper presents a rather different explanation of 
the experimental facts based upon the approximate conservation of baryon- 
minus-lepton number, B — L: If B — L is almost conserved, then the six 
two-component neutrino fields form three nearly Dirac neutrinos, and the 
six neutrino masses coalesce into three nearly degenerate pairs. 

If there are three right-handed neutrinos, then there are six left-handed 
fields, the three left-handed flavor eigenfields u e , and u T and the charge 
conjugates of the three right-handed neutrinos. The neutrino mass matrix 
is then a 6 x 6 complex symmetric matrix A4 which admits a singular- value 
decomposition M. = UMV^ . The singular values are the six neutrino masses 
rrij, and the unitary matrix V' describes the neutrino mixings. 

An angle V is introduced that describes the kind of the neutrinos. Dirac 
neutrinos have 9 V = 0, and Majorana neutrinos have 6 U = ir/2. If all Majo- 
rana mass terms vanish, that is if 9 V = 0, then the standard model conserves 
B — L, which is a global U(l) symmetry. It is therefore natural in the sense 
of 't Hooft [2] to assume that 6 V ~ so that this symmetry is only slightly 
broken. The neutrinos then are nearly Dirac fermions and their masses co- 
alesce into three pairs of almost degenerate masses. Thus the approximate 
conservation of B — L explains the tiny mass differences seen in the solar and 
atmospheric neutrino experiments without requiring the neutrino masses to 
be absurdly small. If one sets 6 V ~ 0.003, suppresses inter-generational mix- 
ing, and imposes a quark-like mass hierarchy, then one may fit the essential 
features of the solar, reactor, and atmospheric neutrino data with otherwise 
random mass matrices M. in the eV range. Thus neutrinos easily can have 
masses that saturate the cosmological bound of about 8 eV. Moreover be- 
cause neutrinos are almost Dirac fermions, neutrinoless double-beta decay is 
suppressed by an extra factor ~ sin 2 9 V sin 2 <p u < 10~ 5 , where (j) u is a second 
neutrino angle, and is very slow, with lifetimes in excess of 2 x 10 27 years. 

This B — L model of neutrino masses and mixings leads to these predic- 
tions about future experiments: The three flavor neutrinos oscillate mainly 
into the conjugates of the right-handed fields, which are sterile. Thus all 
experiments that look for the appearance of neutrinos will yield small or null 
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signals, like those of LSND and KARMEN. Secondly because neutrino masses 
are not required to be nearly as small as the solar and atmospheric mass dif- 
ferences might suggest, neutrinos may well be an important part of hot dark 
matter. Thirdly if a suitable experiment can be designed, it should be pos- 
sible to see the tau neutrino disappear. Fourthly, the rate of neutrinoless 
double-beta decay is suppressed by an extra factor ~ sin 2 6 U sin 2 (fi u < 1CT 5 
and hence will not be seen in the Heidelberg/Moscow, IGEX, GENIUS, or 
CUORE experiments. 



Masses and Mixings 

Because left- and right-handed fields transform differently under Lorentz 
boosts, they cannot mix. It is therefore convenient to write the action exclu- 
sively in terms of two-component, left-handed fields. The two-component, 
left-handed neutrino flavor eigenfields v^v^v T will be denoted z/j, for % = 
e,fi,T. The two-component, left-handed fields that are the charge conju- 
gates of the putative right-handed neutrino fields n re , n rfM , n rT will be denoted 
rii = —io 2 njj for % — e, r, where a 2 is the second Pauli spin matrix. 

The six left-handed neutrino fields Vi,rii for % = 1,2,3 can have three 
kinds of mass terms: The fields and rij can form the Dirac mass terms 
iDijViCr 2 nj — iD*jTijCr 2 v}; in a minimal extension of the standard model, the 
complex numbers Dij are proportional to the mean value in the vacuum of 
the neutral component of the Higgs field. The fields n» and rij can form the 
Majorana mass terms iEijnia 2 rij — iE^-nJo^nf, which break B — L. Because 
these mass terms connect right-handed neutrino fields, which are sterile, they 
do not affect neutrinoless double-beta decay, at least in leading order. Within 
the standard model, the complex numbers Eij are simply numbers; in a more 
unified theory, they might be proportional to the mean values in the vacuum 
of neutral components of Higgs bosons. The fields Vi and Vj can form the 
Majorana mass terms iFiju^Uj — iF*jUja 2 uj, which break SU(2)l <S> U(1)y 
and B — L. Because these mass terms connect left-handed neutrino fields, 
they potentially drive neutrinoless double-beta decay. In a minimal extension 
of the standard model, the complex numbers Fij might be proportional to 
the mean values in the vacuum of the neutral component of a new Higgs 
triplet hab = h ba . 
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Since a 2 is antisymmetric and since any two fermion fields x an d anti- 
commute, it follows that x " 2 ^ = V ;o " 2 X an d X^ 2 ^ = y^Xi which implies 
that the 3x3 complex matrices E and F are symmetric E T = E and F T = 
F and that iDijiijcr 2 ^ = iDijVi<j 2 nj. Thus if we introduce the 6x6 matrix 

D E ) w 

and the (transposed) six-vector N T = (z/ e , u T , n e , n^, n T ) of left-handed 
neutrino fields, then we may gather the mass terms into the matrix expression 

-N T Ma 2 N - -N^M*a 2 Nl (2) 

Z Z 

The complex symmetric mass matrix M. is not normal unless the positive 
hermitian matrix A4A4^ is real because When 
the mass matrix M. is real, it may be diagonalized by an orthogonal trans- 
formation. In general M. is neither real nor normal; but like every matrix, it 
admits a singular-value decomposition [3] 

M = UMV ] (3) 

in which the 6x6 matrices U and V are both unitary and the 6x6 matrix M 
is diagonal, M = diag(mi, m 2 , m 3 , m 4 , m 5 , m 6 ), with singular values rrij > 0, 
which will turn out to be the masses of the six neutrinos. 

The free, kinetic action density of a two-component left-handed spinor ip 
is iip^ (d — a • V) ip. Thus by including the mass terms (2), one may write 
the free action density of the six left-handed neutrino fields N as 

C = iN^ (d -a-V)N + l -N T M(i 2 N - l -N^M*a 2 N^ (4) 

Z Z 

from which follow the equations of motion for N 

(d - a- V)N = M*a 2 N^ (5) 

and iV" 1 " 

(do + a ■ V) a 2 N^ = -MN. (6) 

Applying (d + a ■ V) to the field equation (5) for N and then using the field 
equation (6) for , we find that 

(□ - M?M) N = 0, (7) 
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in which we used the symmetry of the matrix M to write M* as M^. 

The singular-value decomposition M. = U MW allows us to express this 
equation (7) in the form 

(□ - M 2 ) V^N = 0, (8) 

which shows that the singular values mi of the mass matrix M. are the 
neutrino masses and that the eigenfield of mass rrij is 

"m,=EV#Vi. (9) 

1=1 

The vector N m of mass eigenfields is thus N m = V*N, and so the flavor 
eigenfields N are given by N = VN m . In particular, the three left-handed 
fields Vi for % = e, //, r are linear combinations of the six mass eigenfields, V; L = 
J2%i VijV mj and not simply linear combinations of three mass eigenfields. 



Experimental Constraints 

The four LEP measurements of the invisible partial width of the Z impose [4] 
upon the number of light neutrino types the constraint iVj, = 2.984 ± 0.008. 
The amplitude for the Z production of two neutrinos v m . and v mh to lowest 
order is A(v mj , u nik ) oc Yh=i V* k Vij, and therefore the x-section for that pro- 
cess is a(v mp v mk ) oc | Yh=i VikYij] 2 - The LEP measurement of the number 
N u of light neutrino species thus implies that the sum over the light-mass 
eigenfields is 

E lEOy 2 = 2.984 ±0.008. (10) 

j,k light i=l 

This constraint on the 6x6 unitary matrix V is quite well satisfied if all 
six neutrino masses are light. For in this all-light scenario, the sum is 

6 3 3 3 3 3 

E E E = E E = E 1 = 3 - 2-984 ± 0.008. (11) 

j,k=l i=l i'=l i=l i'=l i=l 
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If the Hubble constant in units of 100 km/sec/Mpc is h ~ 0.65, then the 
conservative upper bound on the neutrino component of hot dark matter, 
fi„ <, 0.2, implies [5, 6] that the sum of the masses of the light, stable two- 
component neutrinos that interact weakly is bounded by 

J2 rrij < 8eV. (12) 

j light 



The lowest-order amplitude for a neutrino Vi to be produced by a charged 
lepton ej, to propagate with energy E a distance L as some light-mass eigen- 
field of mass rrij <C E, and to produce a charged lepton is 

im?. L 

A{y { - Uv ) oc ]T VvjV^T^r. (13) 

i light 

The lowest-order amplitude for the anti-process, Pj — > P^, involves the com- 
plex conjugate of the matrix V 

im?: L 

Ate-P,) oc 2 V^e"^. (14) 

i light 

To lowest order the corresponding probabilities are 

^ (%(ml - m*)L\ 

j j' light \ Z£/ / 

and 

P(p^p.,)oc £ Vy^W' exp f ~' (mj l~ mj)L ) ■ (16) 
j light V ZjC/ / 

If all six neutrinos are light, then in the limit L/E — > these sums are 5^/. 

If for simplicity we stretch the error bars on the Chlorine experiment 
and average over one year, then the solar neutrino experiments, especially 
Gallex and SAGE, see a diminution of electron neutrinos by a factor of about 
one-half: 

Pso\(Ve -> Ve) - ~» ( 17 ) 

which requires a pair of mass eigenstates whose squared masses differ by at 
least ~ 10~ 10 eV 2 [5]. The reactor experiments, Palo Verde and especially 
Chooz, imply that these squared masses differ by less than ~ 10~ 3 eV 2 [5]. 
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The atmospheric neutrino experiments, Soudan II, Kamiokande III, IMB- 
3, and especially SuperKamiokande, see a diminution of muon neutrinos and 
antineutrinos by about one-third: 

2 

PatmK -»■ v v) ~ 3' ( 18 ) 

which requires a pair of mass eigenstates whose squared masses differ by 
fCT 3 eV 2 < \m) -m\\< 10" 2 eV 2 [5]. 



The B-L Model 



When the Majorana mass matrices E and F are both zero, the action density 
(4) is invariant under the U (1) transformation N' = e l9G N in which the 6x6 
block-diagonal matrix G = diag(J, —I) with I the 3x3 identity matrix. The 
kinetic part of (4) is clearly invariant under this transformation. The mass 
terms are invariant only when the anti-commutator 

{M,G} = 2( F Q _° E )=0 (19) 

vanishes. 

This U(l) symmetry is the restriction to the neutrino sector of the sym- 
metry generated by baryon-minus-lepton number, B — L, which is exactly 
conserved in the standard model. A minimally extended standard model 
with right-handed neutrino fields n r i and a Dirac mass matrix D but with 
no Majorana mass matrices, E — F — 0, also conserves B — L. When B — L 
is exactly conserved, i.e., when D^O but E — F — 0, then the six neutrino 
masses mj collapse into three pairs of degenerate masses because the left- 
handed and right-handed fields that form a Dirac neutrino have the same 
mass. 

Suppose this symmetry is slightly broken by the Majorana mass matrices 
E and F . Then for random mass matrices D, E, and F, the six neutrino 
masses rrij will form three pairs of nearly degenerate masses as long as the 
ratio 

9n Tr(E^E + F^F) , . 
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Random Neutrino Mass Matrices 
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Figure 1: The six neutrino masses are plotted against the parameter sva.9 v 
for a set of random 6x6 mass matrices. 

is small. For a generic mass matrix A4, the parameter sin 2 9 V lies between 
the extremes < sin 2 Q v < 1 and characterizes the kind of the neutrinos. 
The parameter sin 2 9 V is zero for purely Dirac neutrinos and unity for purely 
Majorana neutrinos. 

Let us now recall 't Hooft's definition [2] of naturalness: It is natural to 
assume that a parameter is small if the theory becomes more symmetrical 
when the parameter vanishes. In this sense it is natural to assume that the 
parameter sin 2 6 U is small because the minimally extended standard model 
becomes more symmetrical, conserving B — L, when sin 2 6 V = 0. 

In Fig. 1 the six neutrino masses rrij are plotted for a set of mass matrices 
Ai that differ only in the parameter sin^. Apart from sin9 u , every other 
parameter of the mass matrices M. is a complex number z = x + iy in which 
x and y were chosen randomly and uniformly on the interval [— leV, leV]. 
It is clear in the figure that when sin 2 9 V 0, the six neutrino masses rrij 
coalesce into three nearly degenerate pairs. Although the six masses of the 
neutrinos are in the eV range, they form three pairs with very tiny mass 
differences when sin 2 6 U rs 0. 
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Thus the very small mass differences required by the solar and atmo- 
spheric experiments are naturally explained by the assumption that the sym- 
metry generated by B — L is broken only slightly by the Majorana mass 
matrices E and F. This same assumption implies that neutrinos are very 
nearly Dirac fermions and hence explains the very stringent upper limits [5] 
on neutrinoless double-beta decay. Because the masses of the six neutrinos 
may lie in the range of a few eV, instead of being squashed down to the meV 
range by the seesaw mechanism, they may contribute to hot dark matter in 
a way that is cosmo logically significant. This B — L model with sin 2 9 V 
is the converse of the seesaw mechanism. 

If sin 2 9 V = 0, then there are three purely Dirac neutrinos, and the mixing 
matrix V is block diagonal V = diag(w*, v) in which the 3x3 unitary matrices 
u and v occur in the singular- value decomposition of the 3x3 matrix D = 
umv\ If these three Dirac neutrinos are also light, then unitarity implies 
that the sum of the normalized probabilities is unity X^ =e -P(Vi — > = 1. 
If sin 2 9 V = 1, then this sum is also unity by unitarity because in this case the 
mixing matrix for the six purely Majorana neutrinos is also block diagonal 
V = diag(t>i? ve)- But if there are six light, nearly Dirac neutrinos, then 
each neutrino flavor z/j will oscillate both into other neutrino flavor eigenfields 
and into sterile neutrino eigenfields. In this case this sum tends to be roughly 
a half Y7i'= e P{. u i ~ * v v) — \ as long as sin 2 9 V is small but not infinitesimal. 
Because of this approximate, empirical sum rule [7] for i = e and //, the 
only way in which the probabilities P S oi(^e — > ^i) and -P a tm(^V — ► ^i) can 
fit the experimental results (17) and (18) is if inter-generational mixing is 
suppressed so that v e oscillates into n e and so that oscillates into n M . 
In other words, random mass matrices Ai, even with sin^ m 0, produce 
probabilities P S oi(^e — ► ^e) and -P a tm(^V — ► ^v) (suitably averaged respectively 
over the Earth's orbit and over the atmosphere) that are too small. The 
probabilities P so \{y e — > u e ) and P a tm(^ do tend to cluster around (|, |) 

as required by the experiments when inter-generational mixing is severely 
repressed, that is if the singly off-diagonal matrix elements of D, E, and 
F are suppressed by 0.05 and the doubly off-diagonal matrix elements by 
0.0025. 

It is possible to relax the factors that suppress inter-generational mixing 
to 0.2 and 0.04 and improve the agreement with the experimental constraints 
(17) and (18) (while satisfying the CHOOZ constraint) provided that one also 
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requires that there be a quark-like mass hierarchy. The points in Fig. 2 were 
generated by random mass matrices M. with sin 9 V = 0.003 by using CKM- 
suppression factors of 0.2 and 0.04 and by scaling the i, j-th elements of the 
mass matrices E,F, and D by the factor f(i) * f(j) where / = (0.2,1,2). 
Thus the mass matrix M. has the r, r elements that are larger than its /i, \i 
elements and //, \i elements that in turn are larger than its e, e elements. The 
clustering of the probabilities P so i(^ e ~~ * v e) an d -P a tm(^ ~~ * ^) around (|, |) 
in Fig. 2 shows that the experimental results (17) and (18) are satisfied. The 
vector / was tuned so as to nearly saturate the cosmological upper bound 
(12) of about 8eV. 

In this scatter plot, every parameter of each of the 10000 matrices M. is a 
complex number z = x+iy with x and y chosen randomly and uniformly from 
the interval [— leV, leV]. The solar neutrinos are taken to have an energy of 
1 MeV, and the probability (15) is averaged over one revolution of the Earth 
about the Sun. The atmospheric neutrinos are taken to have an energy of 
1 GeV, and the probability (15) is averaged over the atmosphere weighted 
by sec 9* z in the notation of Fisher et al. [5] . The thousands of singular- 
value decompositions were performed by the lapack [3] driver subroutine 

ZGESVD [8]. 

Neutrinoless double beta decay occurs when a right-handed antineutrino 
emitted in one decay n — > p + e~ + V e is absorbed as a left-handed neutrino 
in the another decay u e + n — > p + e~. To lowest order these decays proceed 
via the Majorana mass term — iF* e ula 2 ul. Let us introduce a second angle 
4> v defined by 

sin 2 = Tr(FtF) 

^ Tr(E^E + F^F) V ; 

We have seen that we may fit the experimental data (17) and (18) by as- 
suming that sin^j, ~ 0.003 and by requiring the mass matrices E, F, and D 
to exhibit quark-like mass hierarchies with little inter-generational mixing. 
Under these conditions the rate of Ov/3/3 decay is limited by the factor 

\F ee \ 2 < sin 2 9, wn 2 (/> v ml e , (22) 

in which m Ue is the heavier of the lightest two neutrino masses. Thus the 
rate of Qv(3(3 decay is suppressed by an extra factor ~ sin 2 9 V sin 2 0„ <, 10~ 5 
resulting in lifetimes Ti 0u/3/3 > 2x 10 27 yr. The B—L model therefore explains 
why neutrinoless double-beta decay has not been seen and predicts that the 
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Random Quark-Like Mass Matrices with sinx^ = 0.003 
T 




0.65 



Figure 2: The probabilities P so i(^e — > v e ) and P a tm(^ — > *v) for 10000 
random mass matrices .M all with the parameter sin#j, = 0.003, with inter- 
generational mixing suppressed, and with a quark-like mass hierarchy. 



current and upcoming experiments Heidelberg/Moscow, IGEX, GENIUS, 
and CUORE will not see 0v(3(3 decay. 



Conclusions 



The standard model slightly extended to include right-handed neutrino fields 
exactly conserves B — L if all Majorana mass terms vanish. It is therefore 
natural [2] to assume that the Majorana mass terms are small compared to 
the Dirac mass terms. A parameter sin 2 6 V is introduced that characterizes 
the relative importance of these two kinds of mass terms. When this pa- 
rameter is very small, then the neutrinos are nearly Dirac and only slightly 
Majorana. In this case the six neutrino masses rrij coalesce into three pairs 
of nearly degenerate masses. Thus the very tiny mass differences seen in 
the solar and atmospheric neutrino experiments are simply explained by the 
natural assumption that sin 6^ ~ 0.003 or equivalently that B — L is al- 
most conserved. In these experiments the probabilities P so \(v e — > v e ) an d 
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-fatm^ —* v v) are respectively approximately one half and two thirds. One 
may fit these probabilities with random mass matrices in the eV range by 
requiring the neutrino mass matrices E, F, and D to exhibit quark-like mass 
hierarchies with little inter-generational mixing. 

This B — L model leads to these predictions: 

1. Because sin 2 9 V ~ and because inter-generational mixing is sup- 
pressed, neutrinos oscillate mainly into sterile neutrinos of the same 
flavor and not into neutrinos of other flavors. Hence rates for the ap- 
pearance of neutrinos, P(z/j — > v^) with i ^ i', are very low as shown 
by LSND and KARMEN. 

2. The assumption that sin 2 6 V is very small naturally explains the very 
small differences of squared masses seen in the solar and atmospheric 
experiments without requiring that the neutrino masses themselves be 
very small. Thus the neutrinos may very well saturate the cosmological 
bound, J2j m j ~ 8 eV. In fact the masses associated with the points of 
Fig. 2 do nearly saturate this bound. Neutrinos thus may well be an 
important part of hot dark matter. 

3. The disappearance of v T should in principle be observable. 

4. In the B — L model, the rate of neutrinoless double-beta decay is sup- 
pressed by an extra factor ~ sin 2 6 U sin 2 (f) u <, 10~ 5 resulting in lifetimes 
greater than 2 x 10 27 yr. Thus the current and upcoming experiments 
Heidelberg/Moscow, IGEX, GENIUS, and CUORE will not see 0u(3(3 
decay. 
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